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Abstract 

We give a necessary and sufficient condition on the cost function so that the map 
solution of Monge's optimal transportation problem is continuous for arbitrary smooth 
positive data. This condition was first introduced by Ma, Trudinger and Wang |17[ l2T] for 
a priori estimates of the corresponding Monge- Ampere equation. It is expressed by a so- 
called cost-sectional curvature being non-negative. We show that when the cost function 
is the squared distance of a Riemannian manifold, the cost-sectional curvature yields the 
sectional curvature. As a consequence, if the manifold does not have non-negative sectional 
curvature everywhere, the optimal transport map can not be continuous for arbitrary 
smooth positive data. The non-negativity of the cost-sectional curvature is shown to be 
equivalent to the connectedness of the contact set between any cost-convex function (the 
proper generalization of a convex function) and any of its supporting functions. When the 
cost-sectional curvature is uniformly positive, we obtain that optimal maps are continuous 
or Holder continuous under quite weak assumptions on the data, compared to what is 
needed in the Euclidean case. This case includes the reflector antenna problem and the 
squared Riemannian distance on the sphere. 
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1 Introduction 



Given A, B two topological spaces, a cost function c : A x B — > K, and //o,A*i two probability 
measures respectively on A and i?, Monge's problem of optimal transportation consists in 
finding among all measurable maps T : A —> B that push forward \x§ onto \i\ (hereafter 
7#/i = /ii) in the sense that 

(1) V£CB Borel = ^(T" 1 ^)), 
a (the ?) map that realizes 

(2) ^{/ ;I ,T H ) W ,^^,}. 

Optimal transportation has undergone a rapid and important development since the pioneering 
work of Brenier, who discovered that when A — B — W 1 and the cost is the distance squared, 
optimal maps for the problem (j2J) are gradients of convex functions [T]. Following this result 
and its subsequent extensions, the theory of optimal transportation has flourished, with gener- 
alizations to other cost functions (HUGS], more general spaces such as Riemannian manifolds 
[THj . applications in many other areas of mathematics such as geometric analysis, functional 
inequalities, fluid mechanics, dynamical systems, and other more concrete applications such as 
irrigation, cosmology. 

When A, B are domains of the Euclidean space M™, or of a Riemannian manifold, a common 
feature to all optimal transportation problems is that optimal maps derive from a (cost-convex) 
potential, which, assuming some smoothness, is in turn solution to a fully non-linear elliptic 
PDE: the Monge-Ampere equation. In all cases, the Monge-Ampere equation arising from an 
optimal transportation problem reads in local coordinates 

(3) det(£> 2 - A{x, V0)) = f(x, V0), 

where (x,p) — > A(x,p) is a symmetric matrix valued function, that depends on the cost function 
c(x,y) through the formula 

(4) A(x,p) = —D 2 xx c{x, y) for y such that — V x c(x, y) = p. 

That there is indeed a unique y such that — V x c(x, y) = p will be guaranteed by condition Al 
given hereafter. The optimal map will then be 

x -> y : -V x c(x,y) = V(f)(x). 

In the case ^4 = 0, equation (jH)) was well known and studied before optimal transportation 
since it appears in Minkowsky's problem: find a convex hypersurface with prescribed Gauss 
curvature. In the case of optimal transportation, the boundary condition consists in prescribing 
that the image of the optimal map lies in a certain domain. It is known as the second boundary 
value problem. 
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Until recently, except in the particular case of the so-called reflector antenna, treated by 
Wang [2Zj, the regularity of optimal maps was only known in the case where the cost function 
is the (Euclidean) squared distance c(x,y) = \x — y\ 2 , which is the cost considered by Brenier in 
PQ, for which the matrix A in (JHJ) is the identity (which is trivially equivalent to the case A= 0). 
Those results have involved several authors, among which Caffarelli, Urbas , and Delanoe. An 
important step was made recently by Ma, Trudinger and Wang [TJj, and Trudinger and Wang 
[2Tj . who introduced a condition (named A3 and A3w in their papers) on the cost function 
under which they could show existence of smooth solutions to ©• Let us give right away this 
condition that will play a central role in the present paper. Let A = Q, B = Q' be bounded 
domains of ~R n on which the initial and final measures will be supported. Assume that c belongs 
to C 4 {VL x Q!). For (x,y) G (ft x ft'), (£, v) G W l x W n , we define 

(5) &c(x, y)(£, v) := D 2 pm Aij v k v x (x,p), p = -V x c(x, y). 

Whenever £, v are orthogonal unit vectors, we will say that & c (x, y)(C,, v) defines the cost- 
sectional curvature from x to y in the directions (£, v). Note that this map is in general not 
symmetric, and that it depends on two points x and y. The reason why we use the word sectional 
curvature will be clear in a few lines. We will say that the cost function c has non-negative 
cost-sectional curvature on [fl x fl'), if 

(6) e c (x,y)(£,v)>0 Vf^lGfOxfl'j.Vf^lGrxr, f 1 v. 

A cost function satisfies condition Aw on (fl x fl') if and only if it has non-negative cost-sectional 
curvature on (fl x ft'), i.e. if it satisfies (jEJ)- 

Under condition Aw and natural requirements on the domains Q, Q', Trudinger and Wang 
showed that the solution to (J3J) is globally smooth for smooth positive measures //o, They 
showed that Aw is satisfied by a large class of cost functions, that we will give as examples 
later on. Note that the quadratic cost satisfies assumption Aw. This result is achieved by 
the so-called continuity method, for which a key ingredient is to obtain a priori estimates on 
the second derivatives of the solution. At this stage, condition Aw was used in a crucial way. 
However, even if it was known that not all cost functions can lead to smooth optimal maps, 
it was unclear whether the condition Aw was necessary, or just a technical condition for the 
a-priori estimates to go through. 

In this paper we show that the condition Aw is indeed the necessary and sufficient condition 
for regularity: one can not expect regularity without this condition, and more precisely, if 
& c (x,y)(£,v) < for (x,y) G (ft x fi'),£ _L v G M n , one can immediately build a pair of C°° 
strictly positive measures, supported on sets that satisfy the usual smoothness and convexity 
assumptions, so that the optimal potential is not even C 1 , and the optimal map is therefore 
discontinuous. This result is obtained by analyzing the geometric nature of condition (JHj). Let 
us first recall that the solution of the Monge-Ampere equation is a priori known to be cost- 
convex (in short c-convex), meaning that at each point x G fl, there exist y G ft' and a value 
4> c (y) such that 

-<p\y) - c(x,y) = (j)(x), 
-0 c (y)-c(x , ,y)<0(x'), Wx'eQ. 
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The function —(p c (y) — c(x,y) is called a supporting function, and the function y — > (p c (y) is 
called the cost-transform (in short the c-transform) of 0, also defined by 

(f) c (y) = sup{-c(x,y) - (j)(x)}. 

(These notions will be recalled in greater details hereafter). We prove that the condition Aw 
can be reformulated as a property of cost-convex functions, which we call connectedness of the 
contact set 

(7) For all x G Q, the contact set G^x) := {y : (f) c (y) = —<fi(x) — c(x,y)} 

is connected. 

Assuming a natural condition on Q' (namely its c-convexity, see Definition I2.9J) this condition 
involves only the cost function since it must hold for any C defined through a c-transform. 

A case of special interest for applications is the generalization of Brenier's cost ||x — y\ 2 to 
Riemannian manifolds, namely c(x,y) = ^d 2 (x,y). Existence and uniqueness of optimal maps 
in that case was established by McCann [IH] . and further examined by several authors, with 
many interesting applications in geometric and functional analysis (for example JUGS])- The 
optimal map takes the form x —>■ exp x (V<j)(x)) for a c-convex potential and is called a gradient 
map. Then, a natural question is the interpretation of condition Aw and of the cost-sectional 
curvature in this context. We show that we have the identity: 

Cost-sectional curvature from x to x = 8 • Riemannian sectional curvature at x. 

(We mean there that the equality holds for every 2-plane.) As a direct consequence of the 
previous result, the optimal (gradient) map will not be continuous for arbitrary smooth positive 
data if the manifold does not have non-negative sectional curvature everywhere. Although the 
techniques are totally different, it is interesting to notice that in recent works, Lott & Villani 
[23], and Sturm [20] have recovered the Ricci curvature through a property of optimal transport 
maps (namely through the displacement convexity of some functionals). Here, we somehow 
recover the sectional curvature through the continuity of optimal maps. 

We next investigate the continuity of optimal maps under the stronger condition of uniformly 
positive sectional curvature, or condition As: 

(8) 3C >0: 6 c (x,y,Z,v)>C \Z\ 2 \v\ 2 , V(x, y) G (fl x Q'), G W 1 x W 1 , £ ± v. 

We obtain that the (weak) solution of is C l or C 1,a under quite mild assumptions on the 
measures. Namely, for B r (x) the ball of radius r and center x, fix being bounded away from 0, we 
need (io(B r (x)) = o(r n_1 ) to show that the solution of © is C 1 and fi (B r (x)) = 0(r n ~ p ),p < 1 
to show that it is C 1,Q , for a = a(n,p) G (0, 1). Those conditions allow fio, Hi to be singular 
with respect to the Lebesgue measure and no to vanish. 

This results extends the C 1,a estimate of Caffarelli 0] to a large class of cost functions 
and related Monge-Ampere equations. It also shows that the partial regularity results are 
better under As than under Aw, since Caffarelli's C l,a regularity result required /x , A*i to have 
densities bounded away from and infinity, and it is known to be close to optimal [25] . 
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Finally, we prove that the quadratic cost on the sphere has uniformly positive cost-sectional 
curvature, i.e. satisfies As. We obtain therefore regularity of optimal (gradient) maps under 
adequate conditions. 

The rest of the paper is organized as follows: in section 2 we gather all definitions and 
results that we will need throughout the paper. In section 3 we state our results. Then each 
following section is devoted to the proof of a theorem. The reader knowledgeable to subject 
might skip directly to section 3. 

2 Preliminaries 

2.1 Notation 

Hereafter dVol denotes the Lebesgue measure of K n (or of the unit sphere § n_1 for Theorem 
I3.11J) . and B T (x) denotes a ball of radius r centered at x. For 5 > 0, we set classically Qg = 
{x G Q,d(x,dQ) > 5)}. When we say that a function (resp. a measure) is smooth without 
stating the degree of smoothness, we assume that it is C°°-smooth (resp. has a C°°-smooth 
density with respect to the Lebesgue measure). 

2.2 Kantorovitch duality and c-convex potentials 

In this section, we recall how to obtain the optimal map from a c-convex potential in the 
general case. This allows us to introduce definitions that we will be using throughout the 
paper. References concerning the existence of optimal map by Monge-Kantorovitch duality is 
[lj for the cost \x — y\ 2 , ^2] and JU] for general costs, JH] fo r the Riemannian case, otherwise 
the book j2H] offers a rather complete reference on the topic. 

Monge's problem (j2J) is first relaxed to become a problem of linear programming; one seeks 
now 



where n(/z ,/Ui) is the set of positive measures on M. n x M n whose marginals are respectively 
/io and [ii. Note that the (Kantorovitch) infimum Q is smaller than the (Monge) infimum 
of the cost (0), since whenever a map T pushes forward /i onto /ii, the measure n T (x) : = 
/i (x) <g> S T ( x )(y) belongs to U(fii, Hi)- 

Then, the dual Monge-Kantorovitch problem is to find an optimal pair of potentials (0, ip) 
that realizes 





(10) 



J = sup 




4>(x)dno(x) - / 4>(y)dfii(y); 4>(x) + ij(y) > -c(x,y) 



The constraint on 0, ip leads to the definition of c(c*)-transforms: 
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Definition 2.1 Given a lower semi- continuous function (ft : fl C R™ — > R U {+00}, we define 
its c-transform at y E fl' by 

(11) c (y) = sup{-c(x, y) - 0(x)}. 

Respectively, for f : O' C I" -» K a/so /ou>er semi- continuous, define its c*-transform at x & fl 
by 

(12) tp c *(x) = sup{-c(s, y) - ip(y)}. 

yen' 

A function is said cost-convex, or, in short, c-convex, if it is the c*-transform of another 
function ip : fl' — > R ; i.e. for x G fl, (p(x) = sup yen ,{— c(x, y) — ip{y)}, for some lower 
semi- continuous ip : fl' — > R. Moreover in this case <p cc * := (0 C ) C * = on fl (see l&ffi). 

Our first assumption on c will be: 

AO The cost-function c belongs to C A (fl x fl'). 

We will also always assume that fl, fl' are bounded. These assumptions are not the weakest 
possible for the existence/uniqueness theory. 

Proposition 2.2 Under assumption AO, if fl (resp. fl' ) is bounded, <p c (resp. ip c * ) will be 
locally semi-convex and Lipschitz. 

By Fenchel-Rockafellar's duality theorem, we have I = J . One can then easily show that the 
supremum (fTUj) and the infimum @ are achieved. Since the condition <p(x) + ip(y) > —c(x,y) 
implies ip > <ft c , we can assume that for the optimal pair in J we have ip = <p c and <fi = <ft cc *. 
Writing the equality of the integrals in (0 ITTH) for any optimal 7 and any optimal pair (0, <ft c ) 

we obtain that 7 is supported in + 4> c (y) + c(x, y) — Oj. This leads us to the following 

definition: 

Definition 2.3 (Gradient mapping) Let <fi be a c-convex function. We define the set-valued 
mapping by 

G^x) = {ye fl\<j){x) + <p c {y) = -c(a;,y)}. 

For all x E fl, G^ix) is the contact set between <p c and its supporting function —<fi{x) — c(x, •). 

Noticing that for all y e G^x), 0(-) + c(-,y) has a global minimum at x, we introduce / recall 
the following definitions: 

Definition 2.4 (subdifferential) For cf) a semi-convex function, the subdifferential of (p atx, 
that we denote d(f>(x), is the set 

d(p(x) = jp e R n ,0(y) > (f>(x) + p-(y-x) + o(\x - 
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The sub differential is always a convex set, and is always non empty for a semi-convex function. 

Definition 2.5 (c-subdifferential) If is c-convex, the c- sub- differential of (p at x, that we 
denote d c (p(x), is the set 

d c cj){x) = | - V x c(x,y),y G G^x)}. 
The inclusion 7^ <9 c 0(x) C d<p(x) always holds. 

We introduce now two assumptions on the cost-function, which are the usual assumptions made 
in order to obtain an optimal map. 

Al For all x G f2, the mapping y —>■ — V x c(x, y) is injective on Cl'. Its image is a convex set 
of R n . 

A2 The cost function c satisfies det D\ y c 7^ for all (x, y) G f2 x fV. 
This leads us to the definition of the c- exponential map: 

Definition 2.6 Under assumption Al, for x G £1 we define the c-exponential map at x, which 
we denote by % x , such that 

V(x,y) G (tf x n'),Z x (-V x c(x,y)) = y. 

Remark. The definition c-exponential map is again motivated by the case cost=distance 
squared, where the c-exponential map is the exponential map. Moreover, notice that 

Under assumptions Al, A2, is single valued outside of a set of Hausdorff dimension 
less or equal than n — 1, hence, if /xo does not give mass to sets of Assort dimension less than 
n — 1, G^ will be the optimal map for Monge's problem while the optimal measure in (JSJ) will be 
vr = IIq®8g 4> { x )- So, after having relaxed the constraint that the optimal n should be supported 
on the graph of a map, one still obtains a minimizer that satisfy this constraint. 

Let us mention at this point that Monge's historical cost was equal to the distance itself: 
c(x, y) = \x — y\. One sees immediately that for this cost function, there is not a unique y such 
that — V x c(x, y) = V<p(x), and the dual approach fails. 

We now state a general existence theorem, under assumptions that are clearly not minimal, 
but that will suffice for the scope of this paper, where we deal with regularity issues. 

Theorem 2.7 Let Q, Q' be two bounded domains of~R n . Let c G C A (Q x Q') satisfy assumptions 
A0-A2. Let fio,Hi be two probability measures on Q and Q' . Assume that /j does not give 
mass to sets of Assort dimension less or equal than n — 1. Then there exists a <i/i a.e. unique 
minimizer T of Monge's optimal transportation problem Moreover, there exists c-convex 
on Q such that T = G^ (see \2.3i) . Finally, if ip is c-convex and satisfies G^,#/io = £ti, then 
Vip = V0 dfio a.e. 



S 



2.3 The Riemannian case 



This approach has been extended in a natural way to compact smooth complete Riemannian 
manifolds by Robert McCann in [T8] . We first see that all definitions can be translated into 
the setting of a compact Riemannian manifold, replacing the flat connexion by the Levi-Civita 
connexion of the manifold. The cost is allowed to be Lipschitz, semi-concave. When c(-, •) = 
\d 2 {-, •) with d(-, •) the distance function, the c-exponential map is the exponential map, the 
map will be x — > exp x (V g 4>), the gradient V 9 being relative to the Riemannian metric g. 

2.4 Notion of c-convexity for sets 

Following we introduce here the notions that extend naturally the notions of convexity / 
strict convexity for a set. 

Definition 2.8 (c-segment) Let p —* % x (p) be the mapping defined by assumption Al. The 
point x being held fixed, a c-segment with respect to x is the image by T x of a segment ofW 1 . 

If for t> , V\ G W 1 we have % x {vi) — y^i — 0, 1, the c-segment with respect to x joining y$ to 
2/i will be {ye, 9 G [0, 1]} where ye = % x {9v\ + (1 — 6>)t> ). It will be denoted [yo,yi] x . 

Definition 2.9 (c-convex sets) Let Q,Q' C K n . We say that Q' is c-convex with respect to 
Q if for all yo,y± G f2', x G Q, the c-segment [yo,yi] x is contained in Q'. 

Remark. Note that this can be said in the following way: for all i6 0, the set —V x c(x, Q') 
is convex. 

Definition 2.10 (uniform strict c-convexity of sets) For Q, Q' two subsets ofM. n , we say 
that Q' is uniformly strictly c-convex with respect to Q if the sets {— V x c(x, Q')} x &n are uniformly 
strictly convex, uniformly with respect to x. We say that Q is uniformly strictly c*-convex with 
respect to Q' if the dual assertion holds true. 

Remark. This definition is equivalent, although stated under a different form, with the 
definition given in 

Remarks on the sub-differential and c-sub-differential The question is to know if we 
have for all c-convex on Q, for all x G f2, d<fi(x) = d c <fi(x). Clearly, when <p is c-convex and 
differentiable at x,the equality holds. For p an extremal point of d<p(x), there will be a sequence 
x n converging to x such that <p is differentiable at x n and lim n V0(x„) = p. Hence, extremal 
points of d(p(x) belong to d c (p(x). Then it is not hard to show the 

Proposition 2.11 Assume that Q' is c-convex with respect to Q. The following assertions are 
equivalent: 

1. For all <p c-convex onVt,xE Q, d c <p(x) = d<j)(x). 

2. For all <fi c-convex on Q, x G Q, <9 c 0(x) is convex. 

3. For all <fi c-convex on Q, x G Q, G^x) is c-convex with respect to x. 
4- For all (ft c-convex on Q, x G Q, G^(x) is connected. 
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2.5 The Monge- Ampere equation 

In all cases, for a C 2 smooth c-convex potential such that G^/io = A*i, the conservation of 
mass is expressed in local coordinates by the following Monge- Ampere equation 

(13) det( J DLc(x,G (x)) + J D 2 0) = |det£%,c| 

where pi = d/ij/dVol denotes the density of ^ with respect to the Lebesgue measure. (See [E| 
for a derivation of this equation, or [TT], [13].) Hence, the equation fits into the general form 

©• 

2.6 Generalized solutions 

Definition 2.12 (Generalized solutions) Let : Q — ► R &e a c-convex function. Then 

• is a weak Alexandrov solution to (Q3J) and only if 

(14) /or all Bcil, fi (B) = fn(G^B)). 
This will be denoted by /jq = GY/ix. 

• is a weak Brenier solution to il.ty if and only if 

(15) for all B' C Q', /i^B') = fi (G^(B')). 
This is equivalent to jii = G^^/Iq. 

Alexandrov and Brenier solutions First notice that in the definition (|15|) . fii is deduced 
from fi , while it is the contrary in (fT4*|) . As we have seen, the Kantorovitch procedure (fTUj) 
yields an optimal transport map whenever fi does not give mass to sets of Hausdorff dimension 
less than n — 1. Moreover, the map will satisfy (fi3j) by construction, and hence will be 
a weak Brenier solution to (|13|) . Taking advantage of the c-convexity of one can show that 
whenever [i\ is absolutely continuous with respect to the Lebesgue measure, GJ^i is countably 
additive, and hence is a Radon measure (see fTf\ Lemma 3.4]); then a Brenier solution is an 
Alexandrov solutions. Note that one can consider /i = G*dVol, this will be the Monge- Ampere 
measure of 0. Most importantly, for /i supported in Q, G^#/i = lfydVol does not imply that 
G*dVol = fi Q , except if Q' is c-convex with respect to Q (see [T7j). 

2.7 Cost-sectional curvature and conditions Aw, As 

A central notion in the present paper will be the notion of cost-sectional curvature & c (x, y). 
Definition 2.13 Under assumptions A0-A2, we can define 

(16) &c(xo,y Q )(^u) = D 4 (x,p) -> -c(x,Z X0 (p)) 

1 J x ,y =-V x c(x ,y ) 

When £, v are unit orthogonal vectors, & c (x , yo)(£, v) defines the cost- sectional curvature from 
xq to yo in directions (£, f)- 
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Note that this definition is equivalent to (|5|). 

Remark. The definition ()16|) would allow to define a 4-tensor which could be compared 
to the Riemannian curvature tensor, however its symmetries seem to be inconsistent with the 
symmetries of the Riemannian curvature tensor (the so-called Bianchi identities). Nevertheless, 
all the informations on the curvature tensor are contained in the sectional curvatures, i.e. two 
different curvature tensors can not lead to the same sectional curvature. We are now ready to 
introduce the conditions: 

As The cost-sectional curvature is uniformly positive i.e. there exists Co > such that for 
all (x,y) E (n x Q'), for all (z/,f) e W 1 x R n ) with £ JL v, 

e c (x,y)(Z,is)>C \ti\ 2 \is\ 2 . 

Aw The cost-sectional curvature is non-negative: As is satisfied with Co = 0. 

Remark . Let c*(y,x) : = c(x,y). Writing & c in terms of (x, y) and not (x,p) as it is done 
in jT^j, one checks that if c satisfies Aw (resp. As) then c* satisfies Aw (resp. As with a 
different constant). 



2.8 Previous regularity results for optimal maps 

The regularity of optimal maps follows from the regularity of the c-convex potential solution of 
the Monge- Ampere equation (|1H|) . the former being as smooth as the gradient of the latter. It 
falls thus into the theory of viscosity solutions of fully non-linear elliptic equations [S], however, 
the Monge- Ampere equation is degenerate elliptic. Two types of regularity results are usually 
sought for this type of equations: 

Classical regularity: show that the equation has classical C 2 solutions, provided the mea- 
sures are smooth enough, and assuming some boundary conditions. Due to the log-concavity 
of the Monge- Ampere operator, and using classical elliptic theory (see for instance ^E]), C°° 
regularity of the solution of (j!3|) follows from C 2 a priori estimates. 

Partial regularity: show that a weak solution of ()13j) is C 1 or C 1,a under suitable condi- 
tions. We mention also that W 2,p regularity results can be obtained. 



The Euclidean Monge- Ampere equation and the quadratic cost This corresponds to 
the case where the cost function is the Euclidean distance squared c(x, y) = \x — y\ 2 , for which 
c-convexity means convexity in the usual sense, G^x) = V<p(x), and equation (fTH|) takes the 
following form 

(17) detD 2 ^ Po{x) 



Pi(V0(x))' 

Here again, we have pi = d/^o/dVol, i — 0,1. Classical regularity has been established by 
Caffarelli [21 El El Ej, Delanoe [12] and Urbas [22! ■ The optimal classical regularity result, found 
in [21 Ej, is that for C a smooth positive densities, and uniformly strictly convex domains, the 
solution of (|T7|) is C 2 ' a (fl). Partial regularity results have been obtained by Caffarelli [3111013, 
where it is shown that for /j, , \i\ having densities bounded away from and infinity, the solution 
of (fT?]) is C l,a . Thanks to counterexamples by Wang [23 those results are close to be optimal. 
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The reflector antenna The design of reflector antennas can be formulated as a problem of 
optimal transportation on the unit sphere with cost equal to — log \x — y\. The potential (height 
function) <p : S n_1 — > M + parametrizes the antenna A as follows: A = {x<ft(x), x G S™^ 1 }. Then 
the antenna is admissible if and only if is c-convex on S n_1 for c(x,y) = — log |x — y\, and 
G$(x) yields the direction in which the ray coming in the direction x is reflected. This is the 
first non quadratic cost for which regularity of solutions has been established. Wang [2E1 EZj 
has shown classical C 2 (and hence C°°) regularity of solutions of the associated Monge- Ampere 
equation when the densities are smooth. In a recent work, Caffarelli, Huang and Gutierrez 
[D] have shown C 1 regularity for the solution (i.e. continuity of the optimal map) under the 
condition that the measures po and p\ have densities bounded away from and infinity. 

General costs and the conditions As, Aw Recently an important step was achieved in 
two papers by Ma, Trudinger, and Wang . They gave in the first paper jTJj a sufficient condition 
(As, called A3 in their paper) for C 2 (and subsequently C°°) interior regularity. In the second 
paper [21], they could lower this condition down to Aw (condition A3w in their paper) to 
obtain a sufficient condition for global C 2 (and subsequently C°°) regularity, assuming uniform 
strict c-convexity and smoothness of the domains. Note that the result under Aw recovers the 
results of Urbas and Delanoe for the quadratic cost. 

Theorem 2.14 (H3, [H]) Let f2, fi' be two bounded domains ofW 1 . Let c satisfy A0-A2. Let 
Ho, Hi be two probability measures on fl,Q' having densities po,Pi- Assume that p\ G C 2 (Q) is 
bounded away from 0, p2 G C 2 (Q') is bounded away from 0. 

1. IfQ' is c-convex with respect to Q, and c satisfies As, then for c-convex on Q such that 
GW/io = pi, <f) G C 3 (Q). 

2. IfQ,Q' are strictly uniformly c,c*-convex with respect to each other and c satisfies Aw, 
for (f> c-convex on Q such that G^^po = p\, (ft G C 3 (Q). 

3 Results 

We present some answers to the following four questions: 

1. Is there a sharp necessary and sufficient condition on the cost function which would 
guarantee that when both measures have C°° smooth densities, and their supports satisfy 
usual convexity assumptions, the solution of (fT3|) ( and hence the optimal map) is C°° 
smooth ? 

2. Is there a similar partial regularity result available under a general condition ? 

3. What are the cost-functions for which connectedness of the contact set holds © ? 

4. When the cost is set to be the squared distance of a Riemannian manifold, what is the 
meaning of conditions Aw, As in terms of the Riemannian metric ? 
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3.1 Condition Aw, connectedness of the contact set and regularity 
issues 

Answer to questions 1 and 3: Condition Aw is necessary and sufficient for regularity of optimal 
maps. Moreover Aw is equivalent to the connectedness of the contact set. 

Theorem 3.1 Let c satisfy AO, Al, A2, As, Q' being c-convex with respect to Q. Then, for 
all 4> c-convex on Q, for all x G Q, we have 

d<f>(x) = d c (f){x). 

Moreover, the contact set 

G$(x) = {y : <j>(x) + <j) c (y) = -c(x, y)} 

is c-convex with respect to x, it is equal to % x (d(p). If Q is assumed c*-convex with respect to 
Q! , then 

GAv) = i x ■ <t>( x ) + <f> c (y) = ~ c ( x , y)} 

is c*- convex with respect to y 

Then we have the second theorem under Aw. 

Theorem 3.2 Let c be a cost function that satisfies AO, Al, A2, Q, Q' being smooth uniformly 
strictly c-convex (resp. c*-convex) with respect to each other. The following assertions are 
equivalent. 

1. The cost function c satisfies Aw in Q x Q' . 

2. For /io,/ii smooth strictly positive probability measures in there exists a c-convex 
potential (f> G C^fi) such that G^,#/i = fJ>i- 

3. For /io,/ii smooth satirically positive probability measures in there exists a c-convex 
potential (f> G C oo (0) such that G^^Hq = fx±. 

4- For all (ft c-convex in Q, for all x G Q, d c <p(x) = d(f>(x). 

5. For all c-convex in Q, for all y G Q' , the set {x : <p(x) + (p c (y) = —c(x,y)} is c*-convex 
with respect to y. 

6. C 1 c-convex potential are dense among c-convex potentials for the topology of local uniform 
convergence. 

Hence, if condition Aw is violated at some points (x , y ) G (f2 x Q'), there exist smooth positive 
measures /io, Hi on Q, Q' such that there exists no C 1 c-convex potential satisfying G^#/io = fJ>i- 
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Remark. Setting c*(y,x) = c(x,y) we have seen that (5 C > implies © c * > 0. Hence all 
of those assertions are equivalent to their dual counterpart. 
We can add the following equivalent condition for Aw: 

Theorem 3.3 Under the assumptions of Theorem \3.2\ condition Aw holds if and only if, for 
any xq G fi, (yo,yi) G fi', letting cf) be defined by 

4>(x) = max{-c(x, y ) + c(x , y ), -c(x, y x ) + c(x , y x )}, 

for any ye G [yo,yi]x (see Definition \2.fy) . 

4>(x) > -c(x,y e ) +c(x ,ye) 

holds in Q. 

In other words, fe(x) = —c(x, ye) + c(xo, ye) which the supporting function that interpolates at 
xo (nonlinearly) between fo(x) = —c(x,yo) +c(xo,yo) and f\(x) = —c(x,yi) + c(xo,yi), has to 
remain below max{/ , f\}. 

Remark 1. The function furnishes the counter-example to regularity when Aw is not 
satisfied, since for a suitable choice of x ,y ,yi (ft can not be approximated by C 1 c-convex 
potentials. 

REMARK 2. As shown by Propositions 15. l( 15. 7| a quantitative version of Theorem IH.HI 
holds to express condition As. 

Of course condition Aw looks like a good candidate, since in [21] it is shown shown to be 
a sufficient condition for classical regularity. We show here that it is necessary: if it is violated 
at some point, one can always build a counterexample where the solution to (jlHj) is not C 1 
even with C°° smooth positive measures and good boundary conditions (hence the optimal is 
not continuous). Moreover condition Aw is equivalent to a very natural geometric property of 
c-convex functions. 

3.2 Improved partial regularity under As 

Partial answer to question 2: There is partial (i.e. C 1 and C 1,a ) regularity under As, requiring 
much lower assumptions on the measures than what is needed in the quadratic case. There can 
not be C 1 regularity without Aw. When only Aw is satisfied, we don't know yet, except for the 
case c(x, y) = \x — y\ 2 . 

Let us begin by giving the two integrability conditions that will be used in this result. The 
first reads 

For some p e]n, +oo], C Mo > 0, 

(18) Ho{B e {x)) < C Mo e n(1 ~p } for all e > 0, x G fl 
The second condition reads 

For some / : JR + -> R + with lim/(e) = 0, 

(19) Ho{B £ (x)) < f{e)e n{1 -^ for all e > 0, x G fl 
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In order to appreciate the forthcoming theorem, let us mention a few facts on these integrability 
conditions (the proof of this proposition is given at the end of the paper). 

Proposition 3.4 Let Ho be a probability measure on W 1 . 

1. If Ho satisfies J7%|) for some p > n, Ho satisfies (EJ). 

2. If Ho G L p {VL) for some p > n, Ho satisfies KHfy with the same p. 

3. If n Q G L n (Q), fi Q satisfies Uty . 

4- If Ho satisfies JfTty. Ho does not give mass to set of Hausdorff dimension less or equal than 
n — 1, hence \19\) guarantees the existence of an optimal map. 

5. There are probability measures on Q that satisfy Mty) (and hence hlty ) and that are not 
absolutely continuous with respect to the Lebesgue measure. 

Then our result is 

Theorem 3.5 Let c be a cost function that satisfies assumptions AO, Al, A2, As on (QxQ'), 
Q' being c-convex with respect to Q. Let Ho, Hi be probability measures on fl and fl' , and 4> be 
a c-convex potential on fl such that G^^Ho = Hi- Assume that Hi > mdVol on fl' for some 
m > 

1. Assume that Ho satisfies I7%J) for some p > n. Let a — 1 — ^, (3 = 4n " 2 + a ■ Then for any 
5 > we have 

\\(j)\\c^(n s ) < C, 

and C depends only on 5 > 0, C Mo in HT^) . on the constant Co > in condition As, on 
m and on \\D 2 xy c{-, -^lU^nxfi')' \\D 3 c(-, -)IU°°(fixfi')- 

2. If Ho satisfies (TH^). then (ft belongs to C^-^Qs) and the modulus of continuity of V0 is 
controlled by f in < f7^) . 

As an easy corollary of Theorem IH.5| we can extend the C 1 estimates to the boundary if the 
support of the measure Ho is compactly contained in fl. 

Theorem 3.6 Assume in addition to the assumptions of Theorem VJ.b^ that Ho is supported in 
fl 2 , with fl 2 compactly contained in fl. Then, if Ho satisfies M9\) . <fi G C l {fl 2 ) and if Ho satisfies 
JZ3), <j) G C^iCk), with (3 as m Theorem,\TR 

We mention also the consequence of this result for the reflector antenna problem (see paragraph 
12. 81 for a brief description of the problem). 

Theorem 3.7 let Ho, Hi be probability measures on § ra_1 . Let <ft : § ra_1 — > R + be the height 
function of the admissible antenna that reflects the incoming intensity Ho into the outcoming 
intensity Hi- Then the conclusions of Theorems \S.,% \S.(\ hold for <f>. 
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As in the case of the sphere hereafter, we need an a priori estimate to bound the distance 
d(x, G^(x)) away from 0, otherwise the cost function log \x — y\ becomes singular. The strategy 
for the proof of this Theorem will follow the same lines as the proof of Theorem 13.111 

The integrability conditions on fio, f-i are really mild: we only ask that fii be bounded by 
below, and that Ho(B r ) < r n ~ p for p > 1 (p > 1 yields C 1,a regularity) (see conditions (|T8|) and 
(JIHJ) and the subsequent discussion) . The continuity of the optimal map is also asserted in the 
case /io G L n (that implies which is somehow surprising: indeed D 2 (p G L n does not imply 

G C 1 , but here det(D 2 — A(x, V0)) G L n implies G C 1 . Then, as a consequence with 
obtain Theorem 13.71 that improves the result obtained independently by Caffarelli, Gutierrez 
and Huang jU] on reflector antennas. Moreover our techniques yield quantitative C l,a estimates: 
the exponent a can explicitly computed. Finally, our continuity estimates extends up to the 
boundary (Theorem 13.6)1 . This is achieved through a geometric formulation of condition As. 

A full satisfactory answer would include a general result of partial regularity under condition 
Aw. This result is expected in view of the Euclidean case (since the quadratic cost is really the 
limit case for condition Aw). Note that, in view of counterexamples given in [23], the results 
under Aw can not be as good as under As, and can not be much better than Caffarelli's results 
|Hj that require densities bounded away from and infinity. 



3.3 Conditions Aw, As for the quadratic cost of a Riemannian man- 
ifold 

One sided answer to question 4'- The cost- sectional curvature yields the sectional curvature 

Theorem 3.8 Let M be a C A Riemannian manifold. Let c(x,y) = d 2 (x,y)/2 for all (x,y) G 
M x M. Let & c be given by [fity. Then, for all £, z/ G T X M, 

(5 c (x, x) {y } £) 

8 ■ Sectional Curvature of M at x in the 2-plane (£, v). 



l«-(£-")2 



Hence if Aw (resp, As) is satisfied at (x,x), the sectional curvature of M at x is non-negative 
(resp. strictly positive). 

Corollary 3.9 Let M be a compact Riemannian manifold. If the sectional curvature of M is 
not everywhere non-negative, there are smooth positive measures on M such that the optimal 
map (for the cost function c(x,y) = d 2 (x,y)/2) is not continuous. 

At the end of the proof of Theorem 13.81 we give a counterexample to regularity for a two- 
dimensional manifold with negative sectional curvature. 

This observation closes (with a negative answer) the open problem of the regularity of 
optimal gradient maps when the manifold does not have non-negative sectional curvature ev- 
erywhere. Of course, one wonders whether non-negative sectional curvature implies Aw, and 
I could not answer to this question. It seems that when x ^ y the cost-curvature involves 
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derivatives of the curvature tensor along the geodesic linking x to y. However, there is a partial 
converse assertion in the special case of constant sectional curvature: 

The quadratic cost on the round sphere S n_1 satisfies As 

Theorem 3.10 Let S n_1 be the unit sphere ofW 1 equipped with the round metric g, and Rie- 
mannian distance d. Let c(x,y) = \d 2 {x,y). Then c satisfies As on § n_1 x § n_1 \ {(x, x),x G 
S™- 1 }. 

This Theorem is a corollary of the Proposition 16. 2| which shows that condition As is satisfied 
for any choice of local coordinates. Here the cost function is only locally smooth on S n_1 x 
S n-1 \ {(x, x),x G S n_1 }, and it is not enough to show that As is satisfied: we have to prove 
a priori that G^x) stays away from the cut locus of x, otherwise the Monge-Ampere equation 
might become singular (namely the term |-D xy c| blows up). We use for that a previous a-priori 
estimate obtained with Delanoe [Hj, and we can obtain the 

Theorem 3.11 Let E> n be the unit sphere ofMJ 1 equipped with the round metric g, and Rie- 
mannian distance d, and let c(x,y) = \d 2 {x 1 y). Let (p be a c-convex potential on § ra_1 ; then 
for all x G S 71-1 , d(f>(x) = d c (f>(x), and the set G^x) = {y G S n_1 ,0(x) + (f) c (y) = —c(x,y)} is 
c-convex with respect to x. 

Let fi , /ii be two probability measures on § n_1 . Assume that G^^lq = H\. Assume that 
fii > mdVol, for some m > 0. Then 

1. If fiQ satisfies M9j) . then <p G C 1 (S n_1 ), and as in Theorem \3.5\ the modulus of continuity 
o/V0 depends on f in 41.9)) . 

2. If Ho satisfies / li<5j) for some p > n, then <p G C 1 ' /3 (S n ~ 1 ) with j3 = j3(n,p) as in Theorem 

3. If fio, [ii have positive C 1 ' 1 (resp. C°°) densities with respect to the Lebesgue measure, 
then G C 3 ' ^- 1 ) for every a G [0, l[ (resp. G C 00 ^ 1 ' 1 ).) 

3.4 Examples of costs that satisfy As or Aw 

We repeat the collection of cost that was given in [TJj, and |2T] . 

• c(x,y) = ^/l + \x — y\ 2 satisfies As. 

• c(x,y) = — a/1 — \x — y\ 2 satisfies As. 

• c(x, y) = (1 + \x — y\ 2 ) p l 2 satisfies As for 1 < p < 2, \x — y\ 2 < 

• c(x,y) — \x — y\ 2 + — g{y)\ 2 f, g '■ K n — > K convex (resp. uniformly strictly convex) 
with |V/|, \ Vg\ < 1 satisfies Aw (resp. As). 

• c(x,y) = ±~\x — y\ p ,p 7^ and satisfies Aw for p = ±2 or p = — | (— only) and As for 
-f <p< 1 (- only). 

17 



• c(x, y) = — log \x — y\ satisfies As on M. n x M, n \ {(x, x), x G R n }. 

• As pointed out in [21], the reflector antenna problem (.2.6.) corresponds to the case 
c(x, y) = — log \x — y\ restricted to S n . This cost satisfies As on S n_1 x S n_1 \ {x — y}. 

• As shown in Theorem the squared Riemannian distance on the sphere satisfies As. 
Note that it is the restriction to S n_1 of the cost c(x,y) = 6 2 (x,y), where 9 is the angle 
formed by x and y. 
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4 Proof of Theorems [HZD, EL2 

We begin with the following uniqueness result of independent interest: 

Proposition 4.1 Let fi,i> be two probability measures on with Q and Q' connected do- 

mains of MP. Assume that either fi or v is positive Lebesgue almost everywhere in Q (resp. 
in Q'). Then, among all pairs of functions (4>,ip) such that <fi is c-convex, ip is c*-convex, the 
problem M(^) has at most one minimizer up to an additive constant. 

The proof of this proposition is deferred to the end of the paper. 

4.1 Proof of Theorem 13.1 

We will begin with the following lemma: 

Lemma 4.2 Let be c-convex. Let ((f> t ) e>0 be a sequence of c-convex potentials that converges 
uniformly to 4> on compact sets ofQ. Then, if p = —V x c(xo,y) G d(p(xo), xq G Q, y G Q' , there 
exists a sequence (x e ) e >o that converges to xq, a sequence (y e )e>o that converges to y such that 
p e = —V x c(x t ,y e ) G dcf) € (x e ). Finally, p e converges to p. 

Proof. Let y = 1 XQ (p), i.e. p = —V x c(xo,y). Assuming that fl,Q' are bounded, one can 
choose K large enough so that 

4>(x) := 4>(x) + K\x — xq\ 2 /2 + c(x, y), 
<p e (x) := (f) e (x) + K\x - x \ 2 /2 + c(x, y), 

are convex. One can also assume, by subtracting a constant that <P(xq) = 0, and that (f>(x) > 
on Q. Finally, one can assume (by relabelling the sequence) that on B t {xq) compactly contained 
in Q we have \<f> e — <f>\ < e. 
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Consider then 0f = (f> e + S\x — xo| 2 /2 — e. We have (f> 5 e (xq) < 0, and on dB^xo), with /i < r, 

> 0(2) + V/2 - 2e 

> 5^/2 -2e. 

By taking /1 = e 1//3 ,<5 = 4c 1 / 3 , we get that 0f has a local minimum in B^xq), hence at some 
point x e G B^xo), we have 

d4> £ (x e ) 3 -V x c(x e , y) - K(x e - x Q ) -8(x e - x ). 

Then we have \{K + <5)(x e — xo)| small, and thanks to Al, A2, there exists y e close to y such 
that V x c(x e , y t ) = V x c(x e , y) + K(x t — xq) + 5(x t — xo). Thus — — c(x, y e ) has a critical 
point at x e . This implies that p e = —V x c(x t ,y e ) G d(p e (x € ). Finally, since x t — > x,y e — > we 
conclude p e — > p. 

□ 

Now we prove that <9 c </> = 90. In order to do this, we must show that if <ft is c-convex, if 
—</>(•) — c(-, 2/) has a critical point at Xo, this is a global maximum. 
We first have the following observation: 

Lemma 4.3 Let <p be c-convex. Assume that —<p — c(-,y) has a critical point at xq (i.e. G 
d(j)(xo) + V x c(xo,y)), and that it is not a global maximum. Then <fi is not differentiate at xq. 

Proof. Indeed, —</>(■) — c(-,y) has a critical point at Xo, but we don't have <P(xq) + 
4> c (y) = —c(x ,y). However, there is a point y' such that 4>{xq) + 4> c (y') = —c(x ,y'). Hence, 
{— V x c(x , y), — V x c(x , y')} G d4>(x ), and we have V x c(x ,y) ^ V x c(x ,y') from assumption 
Al. 

□ 

We show the following: 

Proposition 4.4 Assume D holds. Let p = —V x c(xo,y) G d<p(xo) with <fi c-convex. Then 
—</>(•) — c(-,y) reaches a global maximum at xq. 

D C 1 c-convex functions are dense in the set j</> c-convex on f2, G<j,(£i) C f2'| for the 
topology of uniform convergence on compact sets of Cl. 

Proof. Assume the contrary, i.e. that — <f>(xi) — c(xi,y) > — 4>{xq) — c(x ,y) for some 
X\ G fl We use D: there exists a sequence of C 1 c-convex potentials (4> e )e>o that converges to 
4>. We use Lemma H~2l there will exist a sequence (x e ) e>0 such that x e — > x and V(j) e (x e ) — > 
— V x c(xo,y). Let y e be such that V0 e (x e ) = — V x c(x e , y e ). Then y e — > y. Since e is C 1 , by 
Lemma ESI #e> the critical point of — e (-) — c(-,y e ) is necessarily a global maximum. Finally, 
since e converges uniformly to 0, we see that —</>(•) — c(-,y) reaches at xo a global maximum. 

Lemma 4.5 L Lei c satisfy AO-As, /et be c-convex with respect to Q. Then D holds. 
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2. Assume are uniformly strictly c-(c*-) convex with respect to each other, and As is 
replaced by condition Aw. Then D holds. 

Proof. As we will see, this result is implied immediately by [T7| (for point 1.) and [21] 
(for point 2.) combined with Proposition 14.11 Let be c-convex. Denote /ii = G^lndVol. 
Note that from Proposition 14. 1[ is the unique up to a constant c-convex potential such that 
G^lodVol = n\. Consider a sequence of smooth positive densities (/4)e>o in such that 
/ifdVol converges weakly-* to fj,i, and has same total mass than Consider e such that 
G^^lndVol = /x^dVol. From 17] in case 1 and [2T| in case 2, e is C 2 smooth inside Q. Then, 
by Proposition 14.11 up to a normalizing constant, e is converging to 0, and V0 e is converging 
to V0 on the points where is differentiable. □ 

Hence, under the assumptions of Lemma l4.5| d<p{x) = d c <fi(x). In view of Proposition 12.1 1( 
the equality d<p(x) = d c <fr(x) for all 0, x is equivalent to the c-convexity of the set 

= {v ■ 00) + 0%) = -c(x,y)j. 
The proof of Theorem Kill is now complete. 

□ 



4.2 Proof of Theorem 

We now show that if Aw is violated somewhere in (Q x Q'), there will exist a c-convex potential 
for which we don't have <90 = d c <f). In view of Lemma 14.51 and Proposition \4A\ this will imply 
that this potential can not be a limit of C 1 -smooth c-convex potentials. Moreover, using the 
construction done in the proof of Lemma 14.51 this implies that there exists positive densities 
/i ,/ii in Q, fl' such that the c-convex potential satisfying G^#fi = /ii is not C 1 smooth. 

Assume that for some xo G Q, y G Q',p = —V x c(xq, y), for some £, v unit vectors in W 1 with 
£ _L v, one has 



(20) D\ vpv p - Dl fXf c(x,Z x (p)) 



> Na > 0. 



Let y = % xo (p — eu), yi = T xo (p + ev), with e small, and recall that y = 1 X0 {p). Hence y is the 
'middle' of the c-segment [yo, yi] x - Let us define 

(21) 0(x) = max | - c(x,y ) + c(x ,y ), -c(x,y ) +c(x ,?/o)}- 

Remark. This function will be used often in the geometric interpretation of As, Aw. It is the 
"second simplest" c-convex function, as the supremum of two supporting functions. It plays 
the role of (xi, ...,x n ) — > \xi\ in the Euclidean case. 

Note first that £ J_ v implies that £ _L (V x c(xo, yi) — V x c(xo, yo))- Consider near xo a 
smooth curve j(t) such that 7(0) = x , 7(0) = £, and such that for t G [—5, 5], one has 

/o(7(*)) : = -c(l(t),yo) +c(x ,y ) = -c(7(£),yi) +c(x ,y 1 ) =: fi(nf(t)). 
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Such a curve exists by the implicit function theorem, and it is C 2 smooth. On 7, we have 



0=2 (/o + /l) 

since /o = fi on 7. Then we compare \{fo + f\) with — c(x,y) + c(xo,y). By (}20|) we have 
n fe^^o, l/o) + D 2 XfXf c{x , yi)] > Dl fXf c(x , y) + c(e, jV ), 



where c(e,N ) is positive for e small enough. Then of course V x c(xo,y) = |[V x c(xo, yo) + 
V x c(xo, yi)}. Hence we have, for e small enough, 

[-c( 7 (t),y) + c(xo,2/)]-0(7(t)) 
= [- C (7(t), y ) + C (x , y )] - I( /o + /l )( 7 ( t )) 



[£>Lc(^o, Vo) + ^Lc(x , yO] - D 2 xx c(x , y) ■ (j(t) - x ) ■ (j(t) - x )/2 + o(t 2 ) 



D l s xs c ( x o, Vo) + D 2 XfXe c{x , yi) 



Z^c(x ,?/) t 2 /2 + o(t 2 ) 



> c(e,N Q )t 2 /2 + o(t 2 ). 



This will be strictly positive for t £ [—5, 5] \ {0} small enough, and of course the difference 
—(f) — [c(x,y) — c(xo,y)] vanishes at Xq. Obviously, the function (p is c-convex, — </>(•) — c(-,y) 
has a critical point at xo, and this is not a global maximum. Hence D can not hold true, since 
we have seen that D would imply that this situation is impossible. 
The proof of Theorems 13.21 EHU is complete. 

□ . 



5 Proof of Theorem 13.5 
5.1 Sketch of the proof 

The key argument of the proof is the geometrical translation of condition As, and how it will 
imply C l regularity for 0: assume that for c-convex cf>, both —</>(■) — c(-, yo) and —</>(■) — c(-, y±) 
reach a local maximum at x — 0. Hence, — Vc(0,yo) and — Vc(0,?/i) both belong to 8(f)(0). 
From assumption Al, we have Vc(0, y\) 7^ Vc(0, ?/o) ; hence <fi is not differentiable at 0. Consider 
ye the c-segment with respect to x = joining |/ to j/i. Then, as we will see in Proposition 
15.11 condition As implies that the functions {— 4>(x) — c(x,yg)} , 9 £ [e, 1 — e] will also have a 
local maximum at x — 0, and moreover this maximum will be strict in the following sense: we 
will have 

-0(x) - c(x, y e ) < -0(0) - c(0, y e ) - 5\x\ 2 /2 + o(|x| 2 ), 

with 5 > depending on |yi — |/o| and Co > in condition As, and bounded by below for 9 
away from and 1. 
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Then, by estimating all supporting functions to on a small ball centered at 0, we will find 
that G^(I? e (0)) contains a Ce neighborhood of {yg,0 6 [1/4,3/4]}, C > depending on C in 
condition As. This is the Proposition 15.31 

Actually, Proposition 15 . 31 shows that one can find locally supporting functions on some small 
ball around 0. Hence we show that dcj) is big. This is where we use Theorem 13. II to show that 
G*(B e (0)) is large. 

One this is shown, we can contradict the bound on the Jacobian determinant of G^. 
We now enter into the rigorous proof of Theorem 13.51 

5.2 Geometric interpretation of condition As 

The core of the proof of Theorem 13.51 is the following proposition, which a geometrical trans- 
lation of assumption As. Actually, (this will be an important remark for section RHIj) . as we 
will see in Proposition 15.71 the result of this proposition is equivalent to assumption As for a 
smooth cost function. 

Proposition 5.1 For yo,yi G Q' , let (ye)ee[o,i] be the c-segment with respect to xq G joining 
yo to yi, in the sense of Definition ^. Si Let 

4>(x) = max{-c(x, y Q ) + c(x , y ), -c(x, y x ) + c(x , yi)}. 

Then for e e]0, \[, for all 9 & [e, 1 — e], for all \x — x \ < Ce, we have 

4>{x) > -c(x,y e ) + c(x ,y e ) + 5 |2/i - Vo\ 2 d(l ~ 0)\x - x \ 2 - j\x - x \ 3 , 

where So depends on Cq > in assumption As, the bound in assumption A2, 7 depends on 
||c(-, -)||c 3 ; an d C is bounded away from for \y \, \yi\ bounded. 

Proof of Proposition 15.11 Rotating and shifting the coordinate, and subtracting an 
affine function, we can assume that xq = and that 

fo(x) := -c(x, y ) + c(0, y ) = ax 1 - D 2 xx c(0, y ) ■ x ■ x/2 + o(\x\ 2 ), 
:= -c(x, yi ) + c(0, Vl ) = bx l - D 2 xx c(0, y x ) ■ x ■ x/2 + o(|x| 2 ), 
-c(x, y e ) + c(0, y e ) = [6b + (1 - 6)a]x 1 - D 2 xx c(0, y e ) ■ x ■ x/2 + o(|x| 2 ), 

with a < b and where x l is the coordinate of x in the direction e^. This means that yo = 
T (ae 1 ),y 1 = T (6ei). 

Using the general fact that max{/ , fi}> Ofi + (1 — 0)f for < 9 < 1, we have 

4>{x) > {6b + (1 - 6)a)x l 

-(6D 2 xx c(0, y x ) + (1 - 6)D 2 xx c(0, y )) ■ x ■ x/2 + o(\x\ 2 ). 

Then we use the assumption As: 
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Lemma 5.2 Under assumption As, 

-D 2 xx c(0,y e )-x-x < -((l-e)D 2 xx c(O,yo) + 0D 2 xx c(O,yi))-x-x 



5\x\ +A|xi||x|, 



where 



5 = -Co|yi — yo\ 2 9(l — 9),C is given in assumption As , 

A = A |yi -y \ 2 9(l - 0),A depends on ||c(-, -)||c 4 (nxn'), 
and \\[det Dxyc] -1 ]] L oo {nxn')- 

PROOF. We first observe the following property on real functions: Let / : R — > 1R be C 2 , 
with f">a> 0. Then we have, for all t , t x e R, 

Of (t ) + (1 - > f(9t + (1 - 0)fi) + - - t \ 2 . 

We apply this observation to the function 

f:t^-D 2 xx c(x,% x (t ei ))-x'-x' 

where x' is equal to (0, x 2 , .., x n ), and hence x' _L e\. /,From assumption As, this function 
satisfies /" > C |x'| 2 . Then, take to — a, t\ — b (note that yg = % x=0 (9b + (1 — 0)a)), to obtain 
that 

-D 2 xx c(0,y e )-x'-x' < -((1 - e)D 2 xx c(0, y ) + 9D 2 xx c(0, Vl )) ■ x' ■ x' 

--C \x'\ 2 6(l-6)\b-a\ 2 . 

Then note that 
(22) 

is uniformly bounded away from and infinity from assumptions Al, A2. To conclude the 
lemma, we now have to get rid of the terms where x 1 appears. First we note that for a C 2 
function /, we have 

Of {to) + (1 - 9)f(h) - f(9t + (1 - 9)h)\ < l\\f\\ c *0(l - 6)\h - t \ 2 . 
We apply this observation to 

/ : t -»• D lx c ( x , Z x (tei)) -x-x- D 2 xx c(x, Ixite-i)) ■ x' ■ x', 

for which we have \f"\ < 2A |x 1 ||a;|, where A depends on ||c(-, ■) \\ C 4 and D p {p — > % x {p)) = 
[D xy c\~ l . Following then the same lines as above, we conclude the lemma. 



\b — a\ 




bt\ — aei 




\yo - yi 




%=o(bei) - 1 x= 


o(aei)| 
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□ 

Using this lemma, we now have 

(23) 4>{x) > (9b + (l-9)a)x x -Dl x c(0,ye)'X-x/2 

+ 5\x\ 2 - A|si||x| + o(|x| 2 ), 

with 5, A given in Lemma 15.21 We need to eliminate the term — A|xi||x|. In order to do so, 
notice that 

\Dl x c(0,y e ) ■ x ■ x/2 - Dl x c(0,y e/ ) ■ x ■ x/2\ < d\9 - 9'\\x\ 2 , 

where G\ depends on |&— a|, \D xy c\~ x and ||c||c3. Then in (J2BJ), using that A |x| \x x \ < Cq\x\ 2 /A+ 
|x 1 | 2 Aq/C , we obtain 

S\x\ 2 - A\x\\x x \ = \y x -y Q \ 2 9{l-9){C Q \ X \ 2 /2-A Q \x\W\) 

> \y x - y \ 2 9(l - 9)(C \x\ 2 /A - (Ag/Co)|a;T). 

Hence for 9,9' G [0,1], (J23|) can be changed in 

4>{x) > (9'b + (1 - 9')a)x 1 - D 2 xx c(0, y' e ) ■ x ■ x/2 

+5\x\ 2 /2 + {{b - a){9 - 9') - C&^x 1 - C x \9 - 9'\\x\ 2 + o{\x\ 2 ), 

where C 2 = 0(1 — 6)\y x — y \ 2 Al/C . We chose e > small. Taking 9' G [e, 1 — e], we choose 
9 such that 9 = 9' + x 1 C2/(b — a). This is possible if we restrict to < (b — a)C 2 _1 e, (using 
(j22J), note that (b — a)^ 1 is bounded away from for \y x — y \ bounded). Note also that 9 
will depend on x x but not 9'. For the term — C x \9 — #'||x| 2 , we have C\ < C x \b — a\ where C x 
depends on \D xy c\~ x and ||c|| C 3. Thus, G x \9 — 9'\\x\ 2 < C x \b — a\\x\ 2 \x 1 \C 2 /\b — a\ < C x C 2 \x\ 3 . 
Hence we find 

We [e,l-e], 

(f){x) > (9'b + (1 - 9')a)x l - D 2 xx c(0, y<y) -x-x/2 + 5\x\ 2 /2 + o(\x\ 2 ). 

Noticing that all the terms o( |x| 2 ) are in fact bounded by C\D 3 c(-, -)| |x| 3 , and that 9'b + 
(1 — 9')a = —V x c(0,ygi), we conclude the Proposition. 

□ 

5.3 Construction of supporting functions 

We let A^(-B) denote the /^-neighborhood of a set B, and we use the Proposition 15. II to prove 
the following: 

Proposition 5.3 Let be c-convex. There exists C > such that if 

\G4,(x x ) - GV(x )| > C\ Xl - xol 1 / 5 > 0, 
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there exists x m G [xo,a;i], such that, if Af e ([xo, Xt]) C Q, then we have 



K({yo, 0e [«,i-«]}) nn'c G^(B e ( Xm )) 



where 




) 



1/2 



#ere y = G</ } (x ),y 1 = G^xx), {y e }ee[o,i} = [yo,yi) Xm denotes the c- segment from y to y 1 with 
respect to x m . The constant C depends only on k and on the constants in the assumptions 
AO-As ; and is bounded away from and infinity for k bounded away from 0. 

Proof. We assume here that we can find xo and x\ in Q close such that V0(xi) — V0(xo) 
is as large as we want compared to X\ — x . If this does not happen, then <p is C 1 ' 1 . We can 
assume that (f>(x ) = <f>(xx)- Then we note y = G,p(x ),yx = G^(xi). From (JUJ), we have 



with equality at x = xo in the first line, at x — xx in the second line. Using that 4>(xq) = <p(xi), 
the supporting functions —c(x, yo) + c(xq, yo) + 4>(xq) and —c(x, y±) + c(x\, y±) + 4>{x\) will cross 
at some point x m on the segment [jcoi^i]- We might suppose that at this point they are equal 
to 0. Hence 



Wx G tt, -c(x , y ) - (f>(x ) > -c(x, y ) - cf>(x), 
Wx G Q, -c(xi,yx) - 4>(xi) > -c(x, y x ) - 4>(x). 



Hence 



c(x, y ) + c(x Q , y ) + 4>(x ) < 4>(x), 
c(x, yi) + c(xi, yi) + <p(xi) < (f)(x), 



(24) 
(25) 



c{x m , y ) + c(x , y ) + <f>(x ) 
c(x m , yi) + c(x 1 ,y 1 ) + <f>(xi) 







0. 



Lemma 5.4 Under the assumptions made above, we have 



< C 3 \x! - x \\yi - yo 



in the segment [xq,xx\, where C3 depends on ||c(-, •) ||c7 2 (Qxr2') ■ 



Proof. Using ()24I25)) . we have 



H 



4>{xo) < -V x c(x m ,y ) ■ (x -x m ) + ||c||oa|x -x m \ 2 /2 
<P{xi) < -V x c(x m , y x ) ■ (x x - x m ) + ||c|| C 2|a;i - x m \ 2 /2. 
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By semi-convexity, on [io,ii] we have <fi{x) < H + C\x\ — xo| 2 , where C is also controlled by 
||c||c 2 (nxfi') Then we assume that -V x c(x m , y ) ■ (x - x m ) and -V x c(x m ,yi) ■ (x x - x m ) are 
both positive, otherwise we are done, since we have assumed that \x 4 — x$\ 2 is small compared 
to \x\ — xqWui — yo I- This implies that 

-V x c(x m ,y ) ■ {x - x m ) < -V x c(x m ,y ) • Oo - x i) 
-V x c(x m , yi) ■ (a?i - x m ) < -V x c(x m , yi) ■ (x 1 - x ). 

Then we have 

2H < -X7 x c(x m ,y )-(x -x 1 )-V x c(x m ,y 1 )-(x 1 -x ) + \\c\\c2\x -x 1 \ 2 

< \V x c(x m ,y ) - V x c(x m ,yi)||x - xx\ + ||c|| C 2|x - xi\ 2 

< \\c\\ C 2^\xi - x \\yi - 2/0 1 + \xq - a;i| 2 ), 

Recall that we assume that \x\ — xq\ is small compared to 1 2/1 — yo\, otherwise there is nothing 
to prove; this means that \x% — x \ 2 is small compared to \x% — xo\\yi — yo\, and we conclude. 

□ 

We now use Proposition 15.11 ( centered at x m ) that will yield 

<j>(x) > max{-c(x,y ) + c(x m ,yo),-c(x,y 1 ) + c(x m ,y 1 )} 
(26) > -c(x,y e ) + c(x m ,y 9 ) + S 9(l - 9)\y - yi| 2 |x - x m \ 2 - j\x - x m \ 3 . 

for all 9 G [e, 1 — e], \x — x m \ < Ce, and with [ye]ee[o,i] the c-segment with respect to x m joining 
yo to y%. Note that e is small but fixed once for all. 

We want to find supporting functions to on a ball of suitable radius. For that we consider 
a function of the form 

fyi x ) = -c(ar, y) + c(x m , y) + <f>(x m ). 

Of course, this function coincides with at x m . If we have f y <(pon dB r (x m ), then —0 + f y 
will have a local maximum inside B r (x m ), hence, for some point x G B r (x m ), we will have 
— V x c(x,y) G d(j)(x). In view of Theorem 13. 1\ this implies that y G G^(B r (x m )). 
First we have 

-c(x,y) + c(x m ,y) 
= -c(x,y e ) + c(x m ,y e ) 

+ / [V ?/ c(x m , y e + s(y - y e )) - V y c(x, y e + s(y - y 6 ))] ■ (y - ye) ds 

Js=0 

< -c(x,ye) + c(x m ,y e ) + ||^c||£«(n x n')l2/ ~ Vo\\ x ~ x m\- 
We then have 

-c(x, y) + c(x m , y) + (f>{x m ) 

< -c(x,ye) + c(x m ,y e ) + C 4 \y - ye\\x - x m \ + C 3 \xi - x \\yi - y \, 
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where C4 = H-D^ciU^Cfixn'), and we have used Lemma E31 to estimate (p(x m ). Using (|2T)|) . we 
want this to be bounded by 

-c(x,y e ) + c(x m ,y e ) + 5 8(1 - 6)\y - yi| 2 |x - x rn \ 2 - j\x - x m \ 3 

on the set {\x — x m \ = r}, for some r > 0. First we restrict 9 to [1/4, 3/4], (i.e. we take e = 1/4 
in (126(0 . Then we want 

3 

— 5 \y - yi\ 2 r 2 - 7r 3 > C 4 \y - y e \r + C 3 |xi - x \\yi ~ y \- 
lb 

We choose \y — yg\ < C$r\yi — y \ 2 for C 5 small enough (for example C 5 = 5o/(16C 4 )), and the 
above inequality will be satisfied for 

2 1 an \ Xl ~ x °\ 



2/i — 2/o| 



if for this value of r, we have indeed < (<5 /16) |?/o — 2/i 1 2 - If n °t then it means that \y± — 2/0 1 5 < 
2 



4pl — ^0 



and we have therefore a C 1 ' 1 / 5 estimate for (/>. 



Now we assume that this is not the case, and therefore the ratio \ X1 x °] is small. Hence we 
consider a ball of radius 

1/2 

r = C R 





- Xq\\ 


V I2/1 


-Vol J 



centered at x m , where Cq = (I6C3) 1 / 2 : . 

We denote \i = C$r\y\ — y \ 2 . Note that C 5 = <5o/(16C 4 ), hence C 5 is bounded away from 0. 
We denote N^S) the /1 neighborhood of a set S. The functions —c(x,y) + c(x m ,y) + <f>(x m ), 
for y G N^{y e , 6 G [1/4, 3/4]} will be equal to <fi at x m , and will be below <fi on the boundary of 
the ball B r (x m ). This proves Proposition 15.31 

□ 



5.4 Continuity estimates 

Proposition 5.5 Let (p be c- convex with G^Q) C Q' . Assume thatc satisfies AO-As inQxQ', 
and that Q' is c-convex with respect to Q. Then, 

• ifGfdVol satisfies then G ClJP); 

• if CrfdVoL satisfies Hty) . for some p > n, then (p G C^'f(fi), witt P(n,p) as in Theorem 
Remark. We recall that G* is defined in Definition 12. 121 

Proof. Consider Q§ = { x £ fi, <i(x, > 5}. In order to have A/" e ([x , #i]) C f2, it is 
enough to have 

1. x ,xi G Us, 
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2. \x -xi\< 5/2, 

3. e < 5/2. 

Note that, in Proposition I5.3[ either |G</,(xo) — Gt/>{%i)\ < C\xi — Xo] 1 ^ 5 , or e < D| X\ — Xq , 
where D depends on C. Hence, for any 5 > small, for all xq, x\ <G fi<5 such that |xo —x%\ < 5/2 
and |xi — xq\ 2 ^ < 5/(2E) where E depends on C, if \G<j,(xo) — G^{x\)\ > C\x\ — x$\ 1 ^, then 
M e ([xo,Xi]) C Q and Proposition 15.31 applies. We now set 

(27) R s = M{5/2,(5/(2E)f 2 }. 

From Proposition IQ1 we will have N^{y e ,6 E [1/4,3/4]} n Q' C G^{B t {x m )). The volume 
of N fl {y e , 9 e [1/4, 3/4]} n O' is of order 

[yi - yo]^- 1 ~ [yi - - yo| a(n_1) , 

(remember that [2/o,2/i]a: m C f2' by c-convexity of Q' with respect to Q) while B e (x m ) has a 
volume comparable to e n . 

C 1,/3 estimates If the Jacobian determinant of the mapping G^ is bounded, (in other words, 
if G*dVol has a density bounded in L°° with respect to the Lebesgue measure) then the volume 
of G^(B e (x m )) must be of order e n , hence we get that \y% — ?/o| 2n_1 < Ce, where C depends on 
all the previous constants Cj. This implies 

/I IX 1 / 2 

|yi-yol M <C7(f^) , 

\ IS/i — 2/o| / 

hence 

|G>(x x ) - G>(x )| < C\xi - xol 3 ^, 

i 

Cj0 G Cjoc _1 (^)' anc ^ we conc hide, using — W x c(x,yi) = V<p(xi),i = 0, 1, that 

|V0(Xi) - V0(x o )| < - Xol 3 ^ 1 ". 

We can refine the argument: Let F be defined by 

F(V) = sup |Vol(G^(£)), 5 C n a ball of volume v\. 
Let // = G* dVol, we have Vol(G (£)) = hence we have F(V) < 

su Pvoi(s)=v Then, 
by Proposition E3 we have F(\B r (x m )\) > Vol(A^{j/ e ,0 e [1/4,3/4]}), which gives 

(28) F L n C% p^Q ) > Crlx, - x \^\ yi - y f^ 
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with u) n the volume of the n-dimensional unit ball, and C-j depends on all the constants Cj 
above, and it can be checked that Cj is bounded away from 0. Assume that F(V) < CV K for 
some kgI. Note that fi E L p implies the (stronger) bound F(V) = o(U 1_1 / p ), hence we might 
write k = 1 — 1/p for some p g]1, +oo], and the condition 

F{V) < CV 1 - 1 ^ 

is then equivalent to condition (JTSj) for \x. Then we find 

M ~ Vol 2 p < C&\xi - £ | 2 p. 
We see first that we need p > n, then we get, setting a = 1 — n/p, 

Vxi, xo G Sis such that \x\ — x$\ < Rs, 
IG^Xt) - G^(y )\ < C 9 \xi - x \ 4n - 2+a , 

where Rs is given in (|27|1. This yields Holder continuity for G^. Then we use that V0(x) = 
— \7 x c(x, G^(x)) and the smoothness of c to obtain a similar Holder estimate for V0. 

C 1 estimates If we only assume condition (|T9*j) for /i = G^dVol, in other words, if F(V) = 
o( yi-i/n^ we can write F( yj < [/(V^/'*)] 3 " -1 V 1_1/n , for some increasing / : [0, 1] -> R+ with 
limy_ > o/(^') — 0. We then have, as |xi — xo| goes to 0, that goes also to (otherwise 

there is nothing to prove). Using the special form of F in (|28j) . we get 

f n ' X ^io^^)>(Ciii2/i-2/o|) 2n - 1 , 

V li/i-vol/ 

hence we get that — y \ goes to when \x\ — xq\ goes to 0. Let g be the modulus of continuity 
of G^ in Qg, then (7 satisfies 

Vxi, x G fi<5 such that \x% — x \ < Rs, f [ C w ) > Cugiu). 

This yields a uniform control on the modulus of continuity of G^\ Indeed / can be chosen 
increasing, thus invertible, and we write 

u>r\C n g{u)) 9 -^. 

wo 

Note that Cxo was obtained from C5, which is bounded away from 0. Then g(u) < u(u) where 

Z 

u is the inverse of z — * f~ (Cnz)—^— . It is easily checked that lim r _^ + u} { r ) = 0- This shows 

C10 

the continuity of G^. Finally we have V0(x) = — V x c(x, G^x)), and the continuity of V0 is 
asserted. 

□ 

Remark 1. The power (3 = 4n " 2+a is n °t optimal for example if n = l,p — +00, for which 
the C 1 ' 1 regularity is trivial, but note that in order to obtain this bound, we had to assume 
that \yi — 2/ j > \x\ — xq] 1 ^ 5 , and, before, that \x\ — x \ = o(\yi — yo\). Hence the conclusion 
should be: either is C 1 ' 1 , or is C hl/5 or is C l > f) '. Note that (3 < 1/7 for n > 2. 
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Proof of Theorem 13.51 In Theorem 13. 5[ we do not want to assume that \L\ E L 1 (IR n ), hence 
we do not have that CrjjUi = f^o, which would imply fio(B) = ^i{G^{B)). Hence we need the 
following proposition to finish the proof: 

Proposition 5.6 Let be c-convex on Q, with G^{fl) C Q' . Assume that G^^fiQ = 
Assume that fix > mdVol on Q f . Then for all u C Q, we have 

fi-ofa) > mVol(G 4>{ui)) , and hence, G*dVol < — // - 

Proof. In Q' we consider N = {y E Q', 3x x ^ x 2 E fi, G^xx) = G 4> (x 2 ) = y}. Then N = 
{y E fi',0 c is not differentiate at y}. Hence Vol(JV) = 0, and Vol(G^(cu) \ N) = Vol(G (c;)). 
Moreover, G^(G^u) \ N) = G^G^lo) \N)Clu. Hence, 

M") > m(G4>( u ) \ N ) > "iVol(w). 

□ 

Proof of the boundary regularity This part is easy: under the assumptions of Theorem 
13.61 the density /i satisfies (fT8)l with p > n (resp. satisfies dHJ). Hence Theorem 13.51 applies 
and <p E C^(fl) (resp. G C^ c (f))). Since f2 2 is compactly contained in Q, we conclude the 
boundary regularity on Q 2 . This proves Theorem 13.61 □ 
Remark. This proof of the boundary regularity is very simple because we have interior 
regularity even when /xq vanishes. This is not the case for the classical Monge- Ampere equation, 
and the boundary regularity requires that both Q and Q' are convex, and is more complicated 
to establish (see j3]). 

We now show that there is indeed equivalence between assumption As at a point x and the 
conclusion of Proposition 15.11 This is a quantitative version of Theorem 13.31 

Proposition 5.7 Assume that at a point xq for all yo,yx, for yx/ 2 the 'middle' point of [yo,yi] xo , 
we have 

4>{x) > -c(x,yi /2 ) + c(x ,yi/ 2 ) + S \y - yi\ 2 \x - x \ 2 + 0(\x - x \ 3 ) 

with cf) as above. Then the cost function satisfies assumption As at xq with Cq = C5q, for some 
constant C > that depends on the bound in A2. 

Proof. The proof follows the same lines as the proof of Theorem 13.21 and is omitted here. 

□ 

6 Proof of Theorem 13.11 

6.1 Monge-Kantorovitch problem on a Riemannian manifold 

The reader interested in details should refer to ^H]- For M a complete Riemannian manifold, 
compact and without boundary, with distance function d(-,-). From now, we restrict to the 
case c(x,y) = ^d 2 (x,y). Then optimal maps for Monge's problem are gradient maps, i.e. 

T opt {x) = G^ix) := exp x (V g (p(x)), 
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where V s - is the gradient with respect to the Riemannian metric g on M, and (j) is some c- 
convex potential obtained by the Monge-Kantorovitch dual approach. With intrinsic notations, 
it has been established ^3] that the Monge- Ampere equation reads 

(29) Jac (expj (V#) d* [ H (cfe .) + # W) = ^J^^y 

where 7i is the Hessian endomorphism operator. Note that for any point x, and p = G^x), 
q — > c(q,p) + 0(g) has a critical point at g = x, therefore we have in local coordinates 

H(c(p,-) + <P)=g i W ij (c(p,-) + ^. 

We will see hereafter the particular form that Dl x c(x, y) takes in the case of the round 
sphere. Finally, in the Riemannian framework, the definitions of sub-differential and c-sub- 
differential take the following form: 



\9J f ' 



d(f>(x) = |p G % x (M),(j)(exp x (v)) >(j)(x) + (v,p) g + o(\v\ 
8 c (j)(x) = |p G % x (M),4>(x) + c (ex P;c (jo)) = -c(x,exp x (p)) 

where (v,p) g denotes the scalar product on % X (M) with respect to the metric g, \v\g = (v,v) g . 



Strategy of the proof Most of the proof is contained in the following points: 

1- Given //o,A*i satisfying the assumptions of Theorem there exists a constant a such 
that d(x , G ^(x)) < tt — o for all x G § n_1 . Hence, we can reduce locally the problem to an 
Euclidean problem, and the distance function does not become singular on that set. 

2- The assumption As is satisfied by the cost function distance squared on the sphere. 
Once this is established, we proceed as follows: 

Given xo G S n_1 we can build around xq a system of geodesic coordinates on the set 
{x, d(x, x ) < R} for R < Ti. From point 1-, for r small enough, the graph {(x, G^x)), x G 
B r (x )} is included in the set B r (x ) x -B 7r -2r(^o) on which the cost function is C°°. From point 
2- and using [T7J, a C 4 smooth solution to (JT3Jl on B r (x ) will enjoy a C 2 a priori estimate at 
Xq. This estimate will depend only on the smoothness of fio, Hi, on r, and r is small but can 
be chosen once for all. Then the method of continuity allows to build smooth solutions for any 
smooth positive densities. 

Then, the Theorem 13.11 follows, in particular that dcj) = d°4>. 

Then Proposition 15.31 holds on B r (x ) x 5 7r _ 2 r(a ; o)- With a straightforward adaption of 
Proposition 15.61 to the sphere, this yields C 1 estimates in, say, B r / 2 (xo), and the Theorem is 
proved. 



6.2 Reduction of the problem to an Euclidean problem 

We denote by B r (resp. B r (x)) a Riemannian ball of radius r (resp. centered at x), §'' 
denotes the unit sphere of M. n . 
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Uniform distance to the cut locus We show that there exists a subset S% of § n_1 x § n_1 on 
which A0-A2 are satisfied, and such that the graph of Gf {(x, G$(x)), x G is contained 

in S%. This subset S% is defined by 



7r — a 



(30) Sl = {(x,y)eS n - 1 xE n -\d(x,y)< 

where a > depends on some condition on /xo,/xi. 

Proposition 6.1 Let /xo,/xi be two probability measures on S n_1 , let <fi be a c-convex potential 
such that G^fiQ = fi\. Assume that there exists m > such that fix > mdVol and that 
Ho satisfies \19\i . Then there exists a > depending on m and , on f in M9\i . such that 
{(x^G^x))^ G S n_1 } C Si, where is defined in l&U\) . 

PROOF: We use [Hj; in that paper, it was shown, for satisfying G</,#/io = /xi, that we 
have \d<p\ < n — e, and e > depends on ||d//i/dVol||x,°° j| (d/xo/dVol) -1 ||^c». Considering <p c and 
G^c that pushes forward /Xi onto /x , we also have a bound such as \d(f>\ < n — e, where e > 
depends on ||d/x /dVol||£oo || (d/ii/dVol)" 1 

Here we slightly extend this bound to the case where /x satisfies (|19jh and /x^ > mdVol. It 
was shown in [TH Lemma 2] that for all Xi, x-i G § n_1 , 



d(x u x 2 ) 

where x[ is the antipodal point to x\. Hence the set E& = {x, d(x, xx) > 5} is sent by G^ in 
B £ (xi) where 

e = In . 

o 

This implies 

Taking 5 > fixed (for example 6 = tt/2), we have 

where D z is the half-sphere centered at z. Hence 

inf \ 1 7r — d<fi(x)\, \tc — d(j) c (x)\ \ 

> inf |r, 3x G § n ~ l , ^(B^x)) > inf /x (A,)} 
:= a. 

In particular, if /xi satisfies and /xo > mdVol, m > 0, we have 

»i(B± r (x)) < (Vol(B 4r (x)))^ l/n 

for r small enough, hence a > and the uniform distance between G^x) and the cut locus of 
x is asserted. Considering C , the same conclusion holds. □ 
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Construction of a local system of coordinates Given x$ G § n_1 , we consider a system of 
geodesic coordinates around Xq, i.e. given a system of orthonormal coordinates at xq and the 
induced system of coordinates on T 3 , (§ n-1 ), we consider the mapping p G T X0 {E> n ~ 1 ) — > exp (p). 
This mapping is a diffeomorphism from Br(0) C M n_1 to Br(xq) as long as R < it. Then, for 
r small enough, we have 



We now take r = a/3 Hence, B r (0) is sent in S 7r _ 2r (0), and the cost function is C°° on 
B r (0) x 5 7r _2r(0)for r small enough. 

In this case, we have Tixip) — ex Pz(p)> an d a l so 



where D v exp x is the derivative with respect to v of t> — > exp x (t>). Assumption Al is trivial, 
since w is indeed uniquely defined by w = exp^fjo). From (}3*Tj) . assumption A2 is true on any 
smooth compact Riemannian manifold, since in this case Jac(t> — > exp x v ) is bounded by above. 

6.3 Verification of assumption As: Proof of Theorem 13.101 

Here we prove the following, which implies Theorem 13.101 

Proposition 6.2 Let x±, ..,x n _i G f2 C MJ 1 ^ 1 be a system of local coordinates on § n_1 , with <p : 
Q — > § n_1 . For any y G S n_1 , x G il 7 letc(x,y) = ^d 2 ((p(x), y). Thenc(-,-) satisfies assumption 
As on Q, with a constant Cq > 0, with Cq bounded away from for {gij(x), g l i(x), 1 < i,j < 
n — l,x G 0} bounded, where [g*-^^] -1 and gfy are the coefficients of the round Riemannian 
metric on § n_1 in the local coordinates x%, ..,x n _i. 

REMARK. It can seem surprising that we do not chose any system of coordinates around 
y. Indeed, y = T x (p) = exp x (p) for p G T x E> n ~ l , and p is written with the induced system of 
coordinates on T X Q. Moreover, condition As is written in terms of x and p only. Indeed, we 
only need to know D 2 xx c{x, exp x (p)) to check condition As. 

Computations in normal coordinates It has been established in jHj that in a system of 
normal coordinates e±, ..e n _i at x with e\ = p/\p\, we have 



This relies on the fact that Dl x c(x, exp x (p)) = [D p (exp x (p))} 1 D x (exp x (p)), and follows by 
computations of Jacobi fields. Hence, we have, in this system of coordinates, 



Si n {{x,y),x G B r (x )} C B r (x ) x B, 




a-2r- 



(31) 



D 2 xy c(x,y) = [D v exp x ] ^exp^ 1 (y)), 




M(x,p) :=D^ B c(a;,X a! (p)) = /-(l 



r cosr 




smr 



33 



where r = \p\, I is the identity matrix of order n — 1, and /' is the identity matrix on 
vect(e 2 , .., e n _i). So, for a given i> e R™ -1 , 



M(x,p) ■ v • v — \v\ — (1 — 



9 . rcosr._ l9 

- " - — rk> 2 , 

smr 



where ITr is the projection of v on p ± . This can be written, using intrinsic notations 



(H(c{-,y))\ y=cxPx{z) v,v) g = \v\ 2 g - (1 - 



\z 


g cos 


z 


9 


sin 


z 


9 



)|n< 

2 



where (•, -) g denotes the Riemannian scalar product on T x (S n ~ ) and \p\ 2 = (p,p) g . 

Let a, (3 G R n_1 , t & M,, u t = a + 1/3 and r = |u t |. We assume that r < n, and ILr is the 
projection of v on uj~. Assumption As is then equivalent to show 

^ [t -> M(x, a + /3t) • u • u] < -C |/3|>| 2 , 



for all a,(3 e T^S™ , r _L /3, which is equivalent to 

rft 2 1/ 



r cos r , 



smr 



)|IW >C |/?|>| 2 , 



for all a,/9 G M n 1 , r _L /3. Without loss of generality, we assume by changing t in £ — 1 that 
ct _L (3. We have 

Ibr = v - \u t \~ 2 {v ■ u t )u t , 

hence, using v _L /3, 

|n t r| 2 = |v| 2 -2(t;-M t ) 2 |M i r 2 + (vu t ) 2 kr 2 

I |2 / \2| 1-2 

= \v | — (r> • a) |u t | . 
Now, assuming |r>| = 1 and (r, a) g = c, we have to evaluate -jmF, where F = G o r, and 



G(r) 
r(t) = \a + 



dt 2 - 

rcosr c 2 % 
sin r r^ 



A computation shows that 



+ 



a 2 {3 2 1 , . Vl c 2 

— — < — 9— — sinr cosrHl -) 

H I sin r r 2 

1 

+- — 

smr 

a ,4 f 2 



2c 2 

sinr — r cosr)^- 



t 2 p 4 



r 2 1 sin 3 r 



(sinr — r cosr)(l -) 



2c 2 



sin 2 r 



(r — sinr cos r)—r- 



smr 



1 , . ,6c 2 . 

smr — r cosrj— - >. 
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Notice that we have for some C > 0, and for all r£ [0, ir] 



(32) 
(33) 



sin r — r cos r > Cr 3 , 
r — sin r cos r > Cr 3 . 



Hence, all terms are positive, except the last one. The sum of the last two lines has for r € [0, n] 
the sign of 

2(r 2 — r sinr cosr) — 3(sin 2 r — r sin r cos r) 
= — (r sin 2r + 3 cos 2r + 4r 2 — 3) . 

One can check (by two successive differentiations) that this quantity is non-negative for r £ 
[0,7r]. Hence we have, using 



^1 



F > C/3 2 



cp 2 



c 2 a 2 + t 2 (3 2 



a 2 c 2 



'1 - 



+ 



a 2 c 2 



Remember that c = a • v, with \v\ = 1, hence lal > |c| and r > |c|, to conclude that 4jrF > 



dt- 



Hence in a normal system of coordinates at x, we see that As is satisfied at x for any y 
such that d(x,y) < it. Through propositions 15.11 and 15. 7[ we see that condition As is satisfied 
in any system of coordinates, with a constant C that depends on the choice of the coordinates 
as in Proposition 16.21 This proves Proposition 16.21 

□ 



6.4 Conclusion of the proof 

C 2 a priori estimates for smooth solutions Let <fi be a C 4 smooth c-convex potential 
such that G<^#/io = /ii, with //o,£*i having positive C 1,1 smooth densities. Consider a system 
of geodesic coordinates around x . As we have seen, using Proposition 16. 1[ in this system of 
coordinates, for r = er/3, we have 

{(x,G^(x)),xe B ro (0)} c £ ro (0) x B n _ 2ro (0) c S 2 a/3 . 

In this system of coordinates, will satisfy on B ro (0) the Monge- Ampere equation (fT^j) where c 
is <i 2 /2 expressed in our local coordinates. Moreover, c satisfies AO- As on B ro (0) x S 7r _ 2r . /3(0). 

Then from 17, Theorem 4.1], we have a bound on D 2 <f)(xo) that depends only on r$ and on 
the bounds on no, fii. Finally, note ro can be chosen once for all once a is known. 
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Smooth solutions In ^3] it was established that given a C°° smooth c-convex potential 0, 
a C°° smooth positive measure /i , and ji\ = G^#/i , the operator 

F : <\> - F(4>) = G^fto 

was locally invertible in C°° around fi\. Hence the existence, for a given pair of C°° smooth 
positive probability measures /i , A*i of a C°° smooth c-convex potential was granted provided 
one could obtain a-priori estimates for the second derivatives. Indeed, this follows from the 
concavity of the equation and and the well known continuity method (see [T6]). 

We conclude the following: for /x , fii having C 1 ' 1 smooth probability densities, there exists a 
(unique up to a constant) c-convex potential £ C 3 ' a for every a £ [0, 1[ such that G^#/io = 
If moreover /x , /ii are C°°, £ C°°. This concludes the third point of Theorem 13. Ill 

C l regularity for weak solutions Having proved the existence of smooth solutions for 
smooth positive densities, we have the first point of Theorem 13.111 Then, using again the 
system of geodesic coordinates, we have G^,(B ro (0)) C B Ro (0) for some Rq < n. We can now 
apply Proposition 15.31 on B ro (0) x .6^(0) (see the remark after). 

Note that V x c(0, B R (0)) = B R (0), which is strictly convex. Hence, c being C 4 on B r (0) x 
B R (0), for r small enough, V x c(x, B R (0)) is convex for all x £ B r (0), and B R (0) is c-convex 
with respect to B r (0). Hence we can now apply Proposition 15.51 on B r (0) x B R (0). 

Then we can adapt with almost no modification Proposition 15.61 as follows: 

Proposition 6.3 Let <ft : S n_1 — > R be c-convex. Assume that G^#/io = H\. Assume that 
fii > mdVol on Q f . Then for all uj C § n_1 , we have 

(i Q (u) > mVol(G^(u;)). 

Hence, we can conclude in the same way the continuity of G^, V0. This achieves the proof of 
Theorem 13.111 

□ 



7 Proof of Theorem I3TT 



For this case it has already been checked (see [21] for example) that the cost function c(x, y) = 
— log |x — y\ satisfies As for x ^ y. We will just prove that under some assumptions on 
the measures //o,£*i we can guarantee that G^x) stays away from x. Then, the proof of the 
Theorem 13.71 follows the same lines as the proof of Theorem 13.111 We prove the 

Proposition 7.1 Let §™ -1 be unit sphere ofW 1 . Let c(x,y) = — \og\x — y\. Let T : S n_1 — > 
§ rt_1 be a 2-monotone map, i.e. such that 

(34) Vxi,x 2 £ S n-1 , c(xi,T(xi)) + c(x 2 ,T(x 2 )) < c(x 2 ,T(x 1 )) + c(x 1 ,T(x 2 )). 

Let n , ii\ be two probability measures on W 1 . Assume that /ii > mdVol for m > 0, and that fj, 
satisfies Utyl . Then there exists eo > depending only on m, f in Uty such that 

Vx £ S n -\d(x,T(x)) > e . 
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Proof. We follow the same lines as in j!4j . From (|34)l . we have 

log \xi - T(x 2 )\ < log \xi - T(xi) \ + log 2 - log \x 2 - T(x 1 )\. 

Letting M = —\og\xi — T(x\)\, the set {x : \x — T(xi)\ > 1} is sent by T in the set {y : 
\y — x i\ < 2 exp(— M)}. Then the bounds on fi , fii yield an upper bound on M as in the proof 
of Proposition 16.11 

□ 



8 Proof of Theorem GHH 



We consider condition Aw at (x ,y = x ). We recall that 

&c(x ,x )(£,v) = -D 2 pvPi/ D 2 HXi [{x,p) c(x,Z X0 (p))}. 

for any u, £ in T XQ M. Let us first take a normal system of coordinates at Xq, so that we will 
compute 

Q = D 2 tt D 2 ss [(x,p) -> c(1 x M)^oM)}. 

Let us write a finite difference version of this operator. We first introduce y_ = % xo {~ hv\ y + = 
{ X Xo (hf), x_ = Ta; (— /i^), x + = X xo (/i£). We use the usual second order difference quotient, for 
example 

Dl^x^^xoip)) = j^(c(x + ,x ) -2c(x ,x Q ) + c(x-,x )). 
(Of course we have c{xq,Xq) = 0.) We will have, as h goes to 0, 



\i,j=+,~ i=+, 



Rearranging the terms, we find that the first line is equal to 

^2 i c ( x -h Vj) ~ cOo, Xi) - c(x , yj)}. 

i,j=+,~ 

Each of the terms inside brackets has a simple geometric interpretation: consider the triangle 
with vertices (x , x iy yj) whose sides are geodesies. This is a square angle triangle. If the metric 
is flat, by Pythagoras Theorem, the term inside the brackets is 0. In the general case, a standard 
computation shows that it is equal to — 2k(xo,^, v)h^ + o(h 4 ) where k(xq,£, v) is the sectional 
curvature at xo in the two-plane generated by £, v. Hence, we get that Q = 8k(x, £, v). 

Now to reach the more general formula of Theorem 13 .81 we use the following expansion of 
the distance that Cedric Villani communicated us: 



37 



Lemma 8.1 Let M be a smooth Riemannian manifold. Let 71,72 be two unit speed geodesies 
that leave point xq G M. Let 6 be the angle between 71 (0) and 72(0) (measured with respect to 
the metric), let k be the sectional curvature of M at x in the 2-plane generated by 7i(0),72(0). 
Then we have 

d 2 (7iW,7 2 (i)) = 2(1 - cos(0))(l - /<cos 2 (#/2))t 2 + 0(t 4 )) 2 . 

Then, we obtain easily, following the same lines as in the case looked above that 

& c (x , x )(f , v) = 8k(x , £, i^)m 2 g W\ 2 g ~ (£, 

where (•, -) g , \ ■ \ g denote respectively the scalar product and the norm with respect to g. This 
proves the Theorem. □ 

8.1 Counterexample to regularity for a manifold with negative cur- 
vature 

Consider the two dimensional surface H = {z = x 2 — y 2 } C M 3 , with its canonical Riemannian 
metric. H has negative sectional curvature around 0. Then for r sufficiently small, Q = 
H R B r (0) is c-convex with respect to itself. Consider the function 

= ma,x{~d 2 /2(X,X ),-d 2 /2(X,X 1 )}, 

where Xo = (0, a, — a 2 ),Xi = (0, —a, —a 2 ). Then, as shown by our proof of Theorem 13 .2\ for 
a small enough, no sequence of C 1 c-convex potentials can converge uniformly to on Q. Let 
fi to be the Lebesgue measure of Q, and \i\ = |(5x + $Xi)- We have G^ # /x = Let 
\i\ G C°°(r2) be a positive mollification of \i\ so that its total mass remains equal to 1, and that 
preserves the symmetries with respect to x = and y — 0. Let </> n be such that G^^fio = [in- 
Then, for n large enough, ra is not differentiable at the origin. Indeed, for symmetry reasons, 
belongs to the subdifferential of (f) n at 0, on the other hand, n converges uniformly to 0, 
and we know from the fact that Aw is violated at that —<f) — c(-, 0) does not reach its global 
maximum on Q at 0. 

9 Appendix 

Proof of Proposition UOJ We prove only the last point, the other points being elementary. 
Consider on M. n a measure locally equal to fiQ = C n ~ x ® /i, where C n ~ l is the n — 1-dimensional 
Lebesgue measure, and // is a probability measure on [0, 1] equal to the derivative of the Devil's 
staircase. Then, /i ^ L 1 . On the other hand, for all [a, b] C [0, 1], fi([a, b]) < \b — a\ a , for some 
a G (0, 1). Then, for x = (x\, .., x n ), 

Vo{B r {x)) < Cr n - X ii\x n -r,x n + r]) < Cr n ~ 1+a = Cr n{1 - 1/p) 

for some p > n. Hence /xq ^ L\ oc and /io satisfies (fT8|) for some p > n. 
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Proof of Proposition I4TTI We know (see |23| chapter 2]) that there exists tt a probability 
measure on K n x MJ 1 , with marginals \i and v, and such that 

(p(x)dfj l (x) + ip(x)du(x) = — J c(x,y)dir(x,y), 

and moreover, there exists a c-convex potential such that 

supp(n) C {(x,G$(x)),x e R n }. 

Let us decompose 7r as 7r = fi^^x, where for d\i almost all j x is a probability measure 

on M. n and ^ x is supported in Gi(x). Hence we have 



d/i(x) 



dj x {y){H x ) + 4>(v) - c i x , y)) 



o. 



This implies that for d\x a.e. x, for d^ x a.e. y, we have y G G^(x). Since for d// a.e. x, we have 
y G Gz(x) d^ x a.s., we deduce that for (i/x a.e. x, (and hence for Lebesgue a.e. x, since /i > 
a.e.), we have Gx(x) fl G^x) ^ 0. This implies that V0 = V0 Lebesgue a.e., and that <t> — 4> 
is constant. This shows that is uniquely defined up to a constant. Now the pair ^ c *,-?/> can 
only improve the infimum (jlUJ) compared to (<f),ip), hence it is also optimal. Hence ip c * is also 
uniquely defined up to a constant. If ip is c*-convex, then ip c * c = ip, and ip is thus uniquely 
defined. 

□ 
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